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Abstract 

We consider the non-trivial zeros of the Riemann ("-function and two classes of L-functions; 
Dirichlet L-functions and those based on level one modular forms. We show that there are an 
infinite number of zeros on the critical line in one-to-one correspondence with the zeros of the 
cosine function, and thus enumerated by an integer n. From this it follows that the ordinate of the 
n-th zero satisfies a transcendental equation that depends only on n. Under weak assumptions, we 
show that the number of solutions of this equation already saturates the counting formula on the 
entire critical strip. We compute numerical solutions of these transcendental equations and also its 
asymptotic limit of large ordinate. The starting point is an explicit formula, yielding an approxi¬ 
mate solution for the ordinates of the zeros in terms of the Lambert lU-function. Our approach is 
a novel and simple method, that takes into account argL, to numerically compute non-trivial zeros 
of L-functions. The method is surprisingly accurate, fast and easy to implement. Employing these 
numerical solutions, in particular for the (-function, we verify that the leading order asymptotic 
expansion is accurate enough to numerically support Montgomery’s and Odlyzko’s pair correla¬ 
tion conjectures, and also to reconstruct the prime number counting function. Furthermore, the 
numerical solutions of the exact transcendental equation can determine the ordinates of the zeros 
to any desired accuracy. We also study in detail Dirichlet L-functions and the L-function for the 
modular form based on the Ramanujan r-function, which is closely related to the bosonic string 
partition function. 
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I. INTRODUCTION 


Riemann’s major contribution to number theory was an explicit formula for the arithmetic 
function 7r(x), which counts the number of primes less than x, in terms of an inhnite sum 
over the non-trivial zeros of the ((s) function, i.e. roots p of the equation ((p) = 0 on the 
critical strip 0 < 3fJ(s) < 1 [1]. It was later proven by Hadamard and de la Vallee Poussin 
that there are no zeros on the line 3fJ(s) = 1, which in turn proved the prime number 


theorem, tt{x) ~ Li(x). (See section VIIC for a review.) Hardy proved that there are 
an inhnite number of zeros on the critical line 3fJ(s) = The Riemann Hypothesis (RH) 
was Riemann’s statement in his seminal eight-page paper [1] that all non-trivial zeros have 
3?(p) = |. In his own words, concerning the roots t to the equation C(| + = 0) 

. .it is very likely that all roots are real. One would of course like to have a rigorous 
proof of this, but I have put aside the search for such a proof after some fleeting vain 

attempts . .. ”. 

Despite strong numerical evidence of its validity, it remains unproven to this day. Many 
important mathematical results were proven assuming the RH, so it is a cornerstone of 
fundamental mathematics. Some excellent introductions to the RH are [2H5]. 


Riemann also gave an estimate N{T), given by (15) but without the S{T) term, for the 


average number of zeros on the entire critical strip with 0 < Q'(p) < T. This formula was 
later proven by von Mangoldt, but it has never been proven to be valid on the critical line, as 
explicitly stated in Edward’s book [2] . Denoting the number of zeros on the critical line up 
to height T by iVo(T), Hardy and Littlewood proved that Nq(T) > CT. Selberg improved 
this result stating that iVo(T) > UTlogT for very small C. Levinson [6] demonstrated that 
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No(T) > CN(T) where C = ^, which was further improved by Conrey [7j who obtained 
C = |. Further improvements on this last result are in Pi- Obviously, if the RH is true 
we must have Nq(T) = N(T). These statements are described in P [TU]. 

The RH is formulated as a problem in pure mathematics, rather than physics, however 
it has interesting connections with different areas of physics such as quantum mechanics, 
quantum chaos, and in particular quantum statistical physics. For an extensive review 
we refer the reader to m- Although the work presented here does not intrinsically bring 
physics ideas to bear on the problem, and is essentially pure mathematics, it is worthwhile 
mentioning some ideas on the RH that are based on physics, even if the purpose is only to 
make contrasts with the present work. Julia and Spector PI proposed independently 
the free “Riemann gas” where the partition function is (. In PI supersymmetry and 
the Witten index were key ingredients. A string theory perspective on the RH is also 
possible. Bakas and Bowick [H] considered an arithmetic gas to construct a formula for 
boson-parafermion equivalence using properties of (. Examples of exactly solvable models 
were also discussed. Spector [15] considered dualities in held theory that are related to 
arithmetic functions. These are analogues of dualities in string theory. He introduced the 
notion of partial supersymmetry, leading to a formulation of parafermions of noninteger order 
and found a bosonic analog of the Witten index. These arithmetic quantum theories have 
a partition function related to (, and possess, like string theory, a Hagedorn temperature. 
In [ini the RH is reformulated in terms of ultraviolet relations occurring in perturbative 
closed strings. A connection between Gromov-Witten invariants, topological string theory 
and Riemann zeros has also been motivated PI- More ideas relating the RH to strings and 
geometry can be found in [18] . A connection of the RH to quantum gases in low dimensions 
was proposed in [T9] . 

The most prominent idea related to physics goes back to an old idea of Hilbert-Polya. 
Below, we will describe and study Montgomery’s conjecture that the ordinates of non-trivial 
zeros of the ^-function satisfy the statistics of the Gaussian Unitary Ensemble (GUE) [20]. 
The latter led Berry to propose that the ordinates of the non-trivial zeros are eigenvalues 
of a chaotic hamiltonian PI- Berry’s work indicates interesting connections of the RH 
to quantum chaos and was further explored in numerous papers. For instance, in [221 [23] 
an analogy between the ordinates of the Riemann zeros and energy levels of a (unknown) 
quantum hermitian operator with chaotic dynamics was proposed. The classical counterpart 
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of such a hypothetical quantum system is associated to the hamiltonian H = xp. Based 
on this approach, a mapping between the Berry-Keating model and the Russian doll model 
of superconductivity was proposed [21]. This model is exactly solvable and has a cyclic 
Renormalization Group. In [25] a generalization of the Berry-Keating model was considered 
by adding an interaction term to the hamiltonian. All these works focus on N{T) and 
carry out the analysis on the critical line, i.e. they essentially assume the validity of the 
RH. Nevertheless, a number of interesting analytic results were obtained, emphasizing the 
important role of the fluctuating term in the counting formula N{T), namely the function 
S{T) = ^ arg( -h iTj. However, these works can only reproduce the smooth part of N{T) 
through a semi-classical approach. Instead of associating Riemann zeros to eigenstates of a 
quantum hamiltonian, as in the previously mentioned papers, the authors of [26] focus on 
the scattering problem. They associate the smooth phase of the ^-function to the density 
of states of a quantum inverted harmonic oscillator. In a related, but essentially different 
approach than Berry and Keating, Connes used abstract mathematical objects called adeles. 
In this approach there exists an operator playing the role of the hamiltonian, which has 
a continuous spectrum, and the Riemann zeros correspond to missing spectral lines IZH. 
Connes proposed a hypothetical trace formula which, if proved, can lead to a prove of the 
RH. A dynamical system whose partition function is the (^-function was also proposed [28] . 
Unfortunately, thus far, a quantum mechanical hermitian operator whose spectrum yields 
the non-trivial zeros has not yet been found. A quantum held theoretical construction with 
a spectrum given by the Riemann zeros has also been pursued, although a free bosonic held 
theory with a spectrum related to prime numbers is unlikely [29l |30] , since its path integral 
cannot be zeta-regularized. We will not be pursuing these ideas here, rather, the basis of 
our work is a novel mathematical analysis of the original problem. 

L-functions are generalizations of the Riemann (^-function, the latter being the trivial 
case [21] • In this paper we will consider two diherent classes of L-functions; Dirichlet L- 
functions and L-functions associated with modular forms. The former have applications 
primarily in multiplicative number theory, whereas the latter in additive number theory. 
These functions can be analytically continued to the entire (upper half) complex plane. 
The Generalized Riemann Hypothesis (GRH) is the conjecture that all non-trivial zeros of 
Dirichlet L-functions and global L-functions in general lie on the critical line. Much less 
is known about the zeros of L-functions in comparison with the ^-function, however let us 
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mention a few works. Selberg [32] obtained the analog of Riemann-von Mangoldt counting 
formula (15) for Dirichlet L-functions. Based on this result, Fuji! [33] gave an estimate 


for the number of zeros in the critical strip with the ordinate between [T,T + H]. The 
distribution of low lying zeros of L-functions near and at the critical line was examined in 
[33], assuming the GRH. The statistics of the zeros, i.e. the analog of the Montgomery- 
Odlyzko conjecture, were studied in [351 - I37] . It is also known that more than half of the 
non-trivial zeros of Dirichlet L-functions are on the critical line ]3H]- For a more detailed 
introduction to L-functions see [33] . 

Besides the Dirichlet L-functions, there are more general constructions of L-functions 
based on arithmetic and geometric objects, like varieties over number helds and modular 
forms [3Ql E]. Some results for general L-functions are still conjectural. For instance, 
it is not even clear if some L-functions can be analytically continued into a meromorphic 
function. We will only consider the additional L-functions based on modular forms here. 
Thus the L-functions considered in this paper have similar properties, namely, they possess 
an Euler product, can be analytic continued into the (upper half) complex plane, except for 
possible poles ai z = 0 and z = 1, and satisfy a non-trivial functional equation. 

Since it is well known that there are an inhnite number of zeros on the critical line for the 
Riemann ^-function, if in some region of the critical strip one can show that the counting 


formula (15) correctly counts the zeros on the critical line, then this proves the RH in this 
region of the strip. It has been shown numerically that the hrst billion or so zeros all lie 
on the critical line [321133] , thus one can approach this problem asymptotically. Such an 
analysis was carried out in [33], where the main outcome was an asymptotic transcendental 
equation for the ordinate of the n-th Riemann zero on the critical line. The way in which 
this equation is derived shows that these zeros are in one-to-one correspondence with the 
zeros of the cosine function; it is in this manner that the n-dependence arises. In this 
paper we provide a more rigorous and through analysis of this result. Moreover, we propose 
generalizations. We derive an exact equation satished by the Riemann zeros on the critical 
line, where the above mentioned asymptotic equation is obtained as a limit of large n. We 
also generalize these results to Dirichlet L-functions and to L-functions related to modular 
forms. For all these classes of functions we obtain an exact equation for the ordinate of 
the n-th zero on the critical line. Since such an equation comes from a relation with the 
cosine function, its solutions can be automatically counted. We will argue that, under weak 
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assumptions, the number of solutions of the transcendental equation coincide with the known 
counting formula for zeros on the entire critical strip, i.e. No{T) = N{T). 

We organize our work as follows. In Section we derive an exact equation satished by 
each individual Riemann zero on the critical line. We discuss how the number of its solutions 


can be the same as the counting formula on the entire critical strip. In Section III we follow 


the same analysis for Dirichlet L-functions, and in Section IV for L-functions based on level 
one modular forms. In Section |V] we derive a useful approximation for the zeros expressed 


explicitly in terms of the Lambert W-function. In Section VI we consider the counterex¬ 
ample of Davenport-Heilbronn, which is known to violate the RH, and discuss how the RH 
fails based on the different properties of our transcendental equation in comparison with 


previous cases. In Section VII we obtain numerical solutions to the transcendental equation 
related to the Riemann (^-function. We show that the leading order asymptotic approxima¬ 
tion is accurate enough to reproduce the GUE statistics and the prime number counting 
function. Furthermore, we show that solutions to the exact transcendental equation yield 
highly accurate results, up to 500 digit accuracy or more if desired. In Section VIII| we solve 
numerically the transcendental equation related to Dirichlet L-functions, considering two 
explicit examples. We also consider numerical solutions for L-functions based on modular 
forms, in particular for the L-function based on the Ramanujan r-function, which is related 
to the bosonic string theory. Section K contains our concluding remarks. In Appendix [X| we 
present the short Mathematica code we used to calculate the zeros for Dirichlet L-functions, 
some of which are shown in Appendix 


II. ZEROS OF THE RIEMANN C-FUNCTION 


For simplicity we hrst consider the Riemann (^-function, which is the simplest Dirichlet 


L-function. Moreover, we hrst consider the asymptotic equation (13), hrst proposed in [H], 
since it involves more familiar functions. However, this asymptotic equation should here be 


viewed as following straightforwardly from the new exact equation (20), presented later. 
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A. Asymptotic equation satisfied by the n-th zero 


Let us start with the completed Riemann zeta function defined by 

xw = x-/n(s/2)c(s) (1) 

where s = a + it. In quantum statistical physics, this function is the free energy of a gas 
of massless bosonic particles in d spatial dimensions when s = d + 1 [19], up to the overall 
power of the temperature Under a “modular” transformation that exchanges one 

spatial coordinate with Euclidean time, if one analytically continues d, physical arguments 
shows that it must have the symmetry 

X(s) = x(l-s). (2) 

This is the fundamental, and amazing, functional equation satisfied by the (^-function, which 
was proven by Riemann using only complex analysis. For several different ways of proving 
(|^ see [T0|. Now consider Stirling’s approximation r(s) ~ which is valid for 

large t. Under this condition we also have 

s® = exp |i ^tlogt + + alogt - ^ + a + O {t~^) | . (3) 

Therefore, using the polar representation ( = and the above expansions, we can 

write x(s) = where 

/ j.\ (o--i)/2 

A(cT,t) = e"^‘/^|C(a + R)| (1+ 0 (s“^)) , (4) 

= \ log - 1) + arg C(cx + it) + O (f^) . (5) 

The above approximation is very accurate. For t as low as 100, it evaluates X (| + R) 
correctly to one part in 10®. Above, we are assuming t > 0. The results for f < 0 follows 
trivially from the relation x(s) = x(s). 

Now let p = cr + R be a Riemann zero. Then arg(^(p) can be defined by the limit 

arg C (p) = lim arg ( {a + 6 + it). (6) 

( 5 —^ 0 + 

For reasons that are explained below, it is important that 0 < 5 1. This limit in general 

is not zero. For instance, for the first Riemann zero given by pi ~ ^ i 14.1347, we have 




argC (pi) 0.157873919880941213041945. On the critical line s = ^ + it, if t does not cor¬ 
respond to the imaginary part of a zero, the well-known function S{t) = ^ arg( -h itJ is 
dehned by continuous variation along the straight lines starting from 2, then up to 2 +it and 
hnally to | -|- it, where argC(2) = 0. Assuming the RH, the current best bound is |*S'(t)| < 
(I -|- o(l)) Jog log ^ for t —)■ cx), proven by Goldston and Gonek pg. On a zero, the more stan¬ 
dard way to dehne this term is through the limit S(p) = | lime_,.o (S (p -|- ie) -|- S' (p — ie)). 
We have checked numerically that for several zeros on the line, our dehnition (|^ gives the 
same answer as this standard approach, and also agrees with the standard dehnition of S(t) 
where t is not the ordinate of a zero. 

From Q we have x(s) = which implies that A{a,—t) = A{a,t) and 6{a,—t) = 

—0{a, t). Denoting y (1 — s) = A' we then have 

A'{a,t) = A{1 — a,t), 6'{a,t) = 6{1 — a,t). (7) 

From (|^ we also have |x(s)| = |x(l —s)|, therefore A{a,t) = A'{a,t) for any s on the critical 
strip. 

Let us now approach a zero p = a + it through the 5 —)■ 0^ limit. From Q it follows that 
({s) and x(s) have the same zeros on the critical strip, so it is enough to consider the zeros 
of x('S)- From ([^ we see that if p is a zero so is 1 — p. Then we clearly have^ 

lim [x(p (5) -h x(l - P - <^)] = lim t)B{a + 5,t) = 0, (8) 

( 5 —>- 0 + 5 —^ 0 + 

where we have dehned 

R(a,t) = + (9) 

The second equality in (|^ follows from A = A'. Then, in the limit 5 —)• O’*", a zero corresponds 
to A = 0, R = 0 or both. They can simultaneously be zero since they are not independent. 
If R = 0 then A = 0, since A oc |C|. However, the converse is not necessarily true. In order 
to be more rigorous, one should consider the limits 5 —)■ O’*" separately in A verses R; below 
we will consider taking the limit in R hrst. 

^ The linear combination in Q was chosen to be manifestly symmetric under s —> 1 — s. Had we taken 
a different linear combination in ([^, then B = for some constant b. Setting the real and 

imaginary parts of B to zero gives the two equations cos0 -I- bcos9' = 0 and sin0 — bsinO' = 0. Summing 
the squares of these equations one obtains cos(0 + 9') = —{b -|- llb)j2. However, since b+l/b > 1, there 
are no solutions except for 5=1. 
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The non-trivial behavior of A is mostly dictated by |C|. On the other hand there is 
much more structure in B since it contains the phases of x(s) and x(l — s). It describes 
oscillations on the complex plane and involves t logt and ( itself. Thus let us consider B = 0. 
We will provide ample evidence that all zeros are characterized by this equation. The general 
solution of 5 = 0 is given by 

e+ e'= {2n + l)7r, (10) 

which are a family of curves t{a). However, since x(s) is analytic on the critical strip, we 
know that the zeros must be isolated points rather than curves, thus this general solution 
must be restricted. Let us choose the particular solution 


6 = 6' lim cos 6 = 0. 

(5—>0+ 


( 11 ) 


On the critical line cr = 1, from ([^ we have that the first equation in ( pT| is already 
satisfied. Then from the second equation in ( [TT| we obtain lim 5 _,.o+ 9 + 6, t) = 

for n = 0, ±1, ±2,..., hence 


n = 


271 


log ( ^ ) -1+ Im - argC (1 + 5 + it) . 

\27ie J 8 <5-5-0+ TT ' 


( 12 ) 


A closer inspection shows that the right hand side of (12) has a minimum in the interval 
(-2,-1), thus n is bounded from below, i.e. n > —1. Establishing the convention that 
zeros are labeled by positive integers, pn = ^ + itn where n = 1 , 2 ,..., we must replace 


n —?• n — 2 in (12). Therefore, the imaginary parts of these zeros satisfy the transcendental 
equation 

t.. f t.. \ 1 11 

(13) 


tyi 1 

2)1 


( + lim 1 argC (1 + i5 + = n - 5^. 

\27re / <5^-0+ TT ' 8 


In short, we have shown that, asymptotically, there are an inhnite number of zeros on the 


critical line whose ordinates can be determined by solving (13). This equation determines 
the zeros on the upper half of the critical line. The zeros on the lower half are symmetrically 
distributed; if = ^ + it^ is a zero, so is = | — it^. 


The left hand side of (13) is a monotonically increasing function of f, and the leading 
term is a smooth function. This is clear since the same terms appear in the staircase 
function N{T), equation (15); see also Remark Possible discontinuities can only come 
from ;^arg(C +it), and in fact, it has a jump discontinuity whenever t corresponds to 
the ordinate of a zero on or off the critical line. However, if lim 5 ^o+ cirg(C (| + 5 + R) is 
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well-defined for every t, then the left hand side of equation (13) is well-dehned for any t, and 
due to its monotonicity, there must be a unique solution for every n. Under this assumption, 


the number of solutions of equation (13), up to height T, is given by 

N„(T) = ^ log (^^) + 1 + i argC (1 + ff) + O (r^). (14) 

This is so because the zeros are already numbered in (13). Thus we can replace n ^ Nq + \ 
and tn —t T, such that the jumps correspond to integer values. In this way T will not 


correspond to the ordinate of a zero and 6 can be eliminated. In summary, Nq(T) in (14) 


counts the solutions to the equation (13) for zeros on the critical line, assuming there is a 
solution for every n, and without assuming the RH. 

Using Cauchy’s argument principle it is known that one can derive the Riemann-von 
Mangoldt formula, which gives the number of zeros inside the critical strip with 0 < %p) < 
T. This formula is given by mm 




(15) 


where S(T) = ;^arg(^ (I +iT). The above formula without the S(T) term was already in 
Riemann’s paper [1]. Note that it has the same form as the counting formula on the critical 


line that we have just found, equation (14). Thus, under the assumptions we have described, 
we conclude that No(T) ~ N(T), at least asymptotically for now. In the next section we 
will present the exact version. This means that our particular solution ( [IT| , leading to 
equation (13), already saturates the counting formula on the entire critical strip and there 


are no additional zeros from A = 0 in nor from the general solution (10). This strongly 


suggests that (13) describes all non-trivial zeros of which must then lie on the critical line. 


We emphasize that we have not assumed the RH in the above arguments. In Section IX 


we will summarize the assumptions which lead to the exact version of the equation (13) 
described below and reiterate its implications for the RH. 
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B. Exact equation for the n-th zero 


It is straightforward to repeat the above analysis without considering an asymptotic 
expansion. The exact versions of Q and (|^ are 

A{a,t) = 7r“'"/2|r (i(a + if)) ||C(a + ff)|, (16) 

6{a,t) = argT + it)) - ^ logvr + argC(cr + if), (17) 

where again x(s) = Ae*® and x(l — s) = with A'{a,t) = A{1 — a,t) and 0'{a,t) = 

9{1 — a,t). The zeros on the critical line o" = ^ correspond to the particular solution 
9 = 9' and lim 5 _,.o+ cos^ = 0. Therefore we have lim 5 _,.o+ 9 + 6,t) = {n+ tt. Replacing 

n —)■ n — 2, the imaginary parts of these zeros must satisfy the exact equation 

argT (I + |f„) -f„log\/^+ lim argC (| + 5 + itn) = (^ - §) (18) 

The Riemann-Siegel d function is dehned by 

d{t) = aigV {\ +It)-t\og^/n, (19) 

where argT is dehned such that this function is continuous and 'd(O) = 0. Therefore, 

we conclude that there are an inhnite number of zeros in the form ^ + itn, where 

n= 1 , 2 ,..., whose imaginary parts exactly satisfy the following equation: 


■9{tn) + lim argC (^ + (5 + itn) = (n - ^) tt. 

< 5 - 5 - 0 + ^ ^ 


( 20 ) 


Expanding the T-function in (19) through Stirling’s formula one obtains 'd(fn) = 
Ylog(^) — I + 0{l/tn), and recovers the asymptotic equation (13) from (20). Let us 


mention at this point that our approach of considering zeros of R = 0, namely (10), is also 
able to reproduce the trivial zeros on the negative real line, and also zeros off of the critical 


line in the counterexample of Section VI [52] . 


Again, as discussed in the paragraph above equation (14), the hrst term in (20) is smooth 
and the whole left hand side is a monotonic increasing function. If lim 5 _,.o+ C (^ + is 


well-dehned for every f, then equation (20) must have a unique solution for every n; see also 


Section VL Under this assumption it is valid to replace tn ^ T and n —)■ Vq + L so the 


number of solutions of (20) is given by 


iVo(T) = -^?(T) + 1 + 1 argC (^ + *T) . 
vr TT 


( 21 ) 
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FIG. 1: (a) A plot of ^ arg^ (5 + it) as a function of t showing its rapid oscillation. The jumps 


occur on a Riemann zero, (b) The function No{T) in (14), which is indistinguishable from a manual 


counting of zeros. The dashed line is the smooth part without the arg term. 

The exact Backlund counting formula, which gives the number of zeros on the entire 
critical strip with 0 < Q'(p) < T, is given by the well-known formula |2] 


N{T) = -i9(T) + 1 + ^ (T). 

TT 


( 22 ) 


Therefore, comparing with the exact counting formula on the entire critical strip (22), we 


have Nq{T) = N{T) exactly. This indicates that our particular solution, leading to equation 


( 20 ), captures all the zeros on the critical strip, and they should all be on the critical line. 


In summary, without assuming the RH, but under the assumption that 
hm 5 _,.o+arg(^ (^-|-(5-f R) exists for every t, then ( 20 ) has a unique solution for every 
n. If one ignores the arg(^ term in the equation, then the unique solution is expressed in 
terms of the Lambert hF-function; see section |Vj If there is indeed a unique solution for 
every n, then this leads to a iVo(T) which saturates the counting formula for the entire 
critical strip, and this would establish the validity of the RH. Furthermore, it implies that 
all non-trivial zeros are simple, as explained in Remark Further related and clarifying 
remarks, based on a counterexample, are in Section |Vl| 


C. Further remarks 


Remark 1. The small shift by 6 in equations (20) or (13) is essential since it smooths out 


S'(t) = b arg ( + it), which is known to jump discontinuously at each Riemann zero. As is 

well known, S{t) is a piecewise continuous function that rapidly oscillates around its average 
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value, which is zero, with discontinuous jumps, as shown in Figure [^. However, when S{t) 
is added to the smooth part of N{T) one obtains an accurate staircase function, which jumps 
by the multiplicity of the zero at the ordinate of each Riemann zero; see Figure [^. Note 
that N{T) is necessarily a monotonically increasing function since it is a counting formula. 


One reason 6 needs to be positive in (20) can be seen as follows. Near a simple 
zero pn we have C(s) ^ (s - Pn) C' (Pn) = {5 + i {t - tn)) C' {Pn)- This gives argC(s) ~ 
arctan ((f — tn)/S) + c, where c is a constant. With 5 > 0 as one passes through a zero from 
below, S(t) increases by one, as it should based on its role in the counting formula N(T). 
On the other hand, if 5 < 0 then S{t) would decrease by one instead, which cannot be the 
case^. 


Remark 2. An important consequence of equation (20) is that, again, if it has a unique 


solution for every n, then all non-trivial zeros are simple. This essentially follows from the 
fact that they are in one-to-one correspondence with the zeros of the cosine function ( [IT| ), 
which are simple. To see this, let us suppose there is a double zero at ordinate t,, i.e. 


tn+i = tn = t,. Then subtracting the equation (20) with tn from the corresponding equation 


with tn+i, we obtain a contradiction, namely 0 = tt. Therefore, for m the equation (20) 
implies tn ^ tm- 

Now if we actually assume that the zeros on the critical line are simple (which we have 
not), there is an easier way to see that the zeros correspond to cos6 = 0. On the critical 
line s = ^ + it, the functional equation ([^ implies x(s) is real, thus for t not the ordinate 
of a zero, sin^ = 0 and cos 6* = ±1. Thus cos 6 is a discontinuous function. Now let t, be 
the ordinate of a simple zero. Then close to such a zero we dehne 

t-t. 


c{t) = 


|x(|+R)| \t-t,\ 


(23) 


For t > t, then c(t) = 1, and for t < t, then c(t) = — 1. Thus c(t) is discontinuous precisely 
at the zero. In the above polar representation, formally c{t) = cos9 Therefore, by 

identifying zeros as the solutions to cos 0 = 0, we are simply dehning the value of the function 
c(t) at the discontinuity as c(t,) = 0. 

Remark 3. Let us introduce another function (C(s) —>■ ^(s) = f{s)C{s) that also satishes the 
functional equation ([^, i.e. x(s) = x(l — s), but has zeros off of the critical line due to the 

^ Note added: there are deeper reasons why 6 has to be positive, described in our subsequent work [KH] . 
which is discussed in the concluding section of the present article. 
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zeros of f{s). In such a case the corresponding functional equation will hold if and only if 
/(s) = /(I — s) for any s, and this is a trivial condition on f{s) which could have been can¬ 
celed in the hrst place. Moreover, if f{s) and C(s) have different zeros, the analog of equation 
@ has a factor /(s), i.e. x{p + S) + x{^-p-^) = /(p + ^) [x(p + 5) + x(l - p - 5)] = 0, im¬ 
plying (|^ again where x(s) is the original Q. Therefore, the previous analysis eliminates 
/(s) automatically and only hnds the zeros of x('S)- The analysis is non-trivial precisely 
because ({s) satishes the functional equation but ({s) ^ C(1 ~ -s)- Furthermore, it is a 
well-known theorem that the only function which satishes the functional equation (|^ and 
has the same characteristics of C(5), is C{s) itself. In other words, if ({s) is required to have 
the same properties of C(-5), then (^(s) = C ({s) where C is a constant [101 pg. 31]. 


Remark 4. Although equations (20) and (22) have an obvious resemblance, it is impossible 
to derive the former from the later, since the later is just a counting formula valid on the 
entire strip, and it is assumed that T is not the ordinate of a zero. Moreover, such a 
derivation would require the assumption of the validity of the RH and the simplicity of the 


zeros, contrary to our approach, where we derived equations (20) and (13) directly on the 
critical line, without assuming the RH, nor the known counting formula N{T). Despite our 


best efforts, we were not able to hnd equations (13) and (20) in the literature. Furthermore, 


the counting formulas (14) and (22) have never been proven to be valid on the critical line 

i- 


Remark 5. One may object that our basic equation (20) involves C(s) itself and this is 
somehow circular. This is not a valid counter-argument. First of all, arg^ already appears 


in the counting function N{T). Secondly, the equation (20) is a much more detailed equation 


than simply C,{s) = 0, which has an inhnite number of solutions, in contrast with (20) which 
for each u, as we have argued, has a unique solution corresponding to the u-th zero. Also, 
there are well known ways to calculate arg for example from an integral representation or 
a convergent series 
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III. ZEROS OF DIRICHLET L-FUNCTIONS 


A. Some properties of Dirichlet Z-functions 


We now consider the generalization of the previous results to Dirichlet L-functions. Let 
us hrst introduce the basic ingredients and dehnitions regarding this class of functions, which 
are all well known [31]. Dirichlet L-series are dehned as 


L{s,x) 


n=l 


(24) 


for 3?(s) > 1, where the arithmetic function x(^) is a Dirichlet character. They can all be 
analytically continued to the entire complex plane, except possibly for a simple pole at s = 1 
when X is principal, and are then referred to as Dirichlet L-functions. 

There are an inhnite number of distinct Dirichlet characters which are primarily char¬ 
acterized by their modulus A;, which determines their periodicity. They can be dehned 
axiomatically, which leads to specihc properties, some of which we now describe. Consider a 
Dirichlet character y mod fc, and let the symbol {n, k) denote the greatest common divisor 
of the two integers n and k. Then y has the following properties: 


1 . xin + k) = xin). 

2 . x(l) = 1 and x( 0 ) = 0 . 

3. xinm) = xin)x{jn). 

4. x(n) = 0 if (n, k) > 1 and x(n) 7 ^ 0 if (n, k) = 1. 

5. If (n. A;) = 1 then = 1, where ‘^{k) is the Euler totient arithmetic function. 

This implies that x(n) are roots of unity. 

6 . If X is a Dirichlet character so is the complex conjugate x- 


For a given modulus k there are ^p{k) distinct Dirichlet characters, which essentially follows 
from Property]^ above. They can thus be labeled as Xkj where j = 1, 2,..., ‘^{k) denotes an 
arbitrary ordering. If A: = 1 we have the trivial character where x(n) = 1 for every n, and 


(24) reduces to the Riemann ^-function. The principal character, usually denoted by Xi-, 
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is defined as Xi(^) = 1 if {n, k) = 1 and zero otherwise. In the above notation the principal 
character is always Xk,i- 

Characters can be classihed as primitive or non-primitive. Consider the Gauss sum 

k 

G{x) = J2x{m)e^^^^/'^. (25) 

m=l 

If the character x mod k is primitive, then |G(x)p = k. This is no longer valid for a non¬ 
primitive character. Consider a non-primitive character x mod k. Then it can be expressed 
in terms of a primitive character of smaller modulus as x(n) = xT(n)y;(n), where ^ is the 
principal character mod k and y is a primitive character mod k < k, where A; is a divisor of 
k. More precisely, k must be the conductor of x (see |3T] for further details). In this case 
the two L-functions are related as L(s, x) = X)^p\k (f ~ x{p)/p^)- Thus L(s, x) h^is the 
same zeros as L{s,x)- Therefore, it suffices to consider primitive characters, and we will 
henceforth do so. 

We will need the functional equation satished by L{s, y). Let y be a primitive character. 
Dehne its order a such that 


if X(-1) = -1 (odd) 
I 0 if y(—1) = 1 (even) 


(26) 


Let us dehne the entire function 

3+a 

(27) 

Then A satishes the following well-known functional equation, only valid for primitive char¬ 
acters pi] : 

A(s,y) = ^—^^A(l - s,y). (28) 


B. Exact equation for the n-th zero 


For a primitive character, since |G(y)| = \/k, the factor on the right hand side of (28) is 


a phase. It is thus possible to obtain a more symmetric form through a new function dehned 
as 

.,-a/2 ^1/4 

(29) 


^(®W) = ^ A(s,y). 
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It then satisfies 


= ^(1 - s,x) = ^(1 - s,x)- (30) 

Above, the function of s is defined as the complex conjugation of all coefficients that define 
namely y and the factor, evaluated at a non-conjugated s. 


Note that A(s,x) = A(s,x). Using the known result G (x) = x(—l)G(x) we then con¬ 
clude that 

= ^(s,x) • (31) 

This implies that if the character is real, when p is a zero of ^ so is p, and one needs only 
to consider p with positive imaginary part. On the other hand if y 7 ^ y, then the zeros 
with negative imaginary part are different than p. For the trivial character where k = 1 


and a = 0, implying y(n) = 1 for any n, then L(s, y) reduces to the Riemann ^(s) and (30) 
yields the well-known functional equation (|^. 

Let s = a + it. Then the function (29) can be written as ^(s, y) = Ae*® where 


-4(a,f,y) = ( - 

TT 


0{cr,t,x) = argT 


(T+g 


a + a + it 


\L{a + it,x)\, 


a + a + it 


'Ml 


2argG'(\') 


(32) 

(33) 


vra 


+ aig L{a + it, x) + ^ 


From (31) we have that A(cr, f,y) = A(cr, —f,y) and 6{a,t,x) = Denoting 

^(1 - s, X) = we then have A' {a, t, y) = A (1 - a, t, y) and 9'{a, t, y) = 9{1 - a, t, y). 


Taking the modulus of (30) we also have that A{a, t, y) = A'{a, t, y) for any s. 

On the critical strip, the functions L{s, y) and ^{s, y) have the same zeros. Thus on a 
zero we clearly have 

lim {^(p + (5,y)+^(l-p-5,y)} = 0 . 

Let us define 


(34) 

(35) 


Since A = A' everywhere, from (34) we conclude that on a zero we have 

lim A(cr -F 6, t, x)B{a + 5, t, y) = 0. 
s^o+ 

As before, let us consider the particular solution of lim 5 _,.o+ B = 0 given by 


(36) 


9 = 9', lim cos 0 = 0. 

(5—>-0+ 


( 37 ) 
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Let us define the function 


n /X -rx /1 a ■ t\ t , 

i9k,ai^) = r log 


= 9 iogr( ^ + ^ + *- 


, t 
2 


'Ml 


(38) 


When k = 1 and a = 0, the function (38) is just the usual Riemann-Siegel "d function 


(19). Since the function logL has a complicated branch cut, one can use the following series 


representation in (38) 153 

OO 

log r(s) = -7S - log s - ^ jlog 


n=l 


n 


1 + - -- 


n 


(39) 


where 7 is the Euler-Mascheroni constant. Nevertheless, most numerical packages already 
have the logL function implemented. 


On the critical line cr = ^ the first equation in (37) is already satished. From the second 


equation we have hm 5 _,.o+ 9 + 5, t) = (n + |) vr, therefore 

n / \ 1 - T n e • \ argG vra / 1 

Vk,a[tn) + ^Im argL (2 + c> + Rn, x) -^-h — = I n + - I TT. 


(40) 


Analyzing the left hand side of (40) we can see that it has a minimum, thus we shift 
n —>■ 77 , — (hq + 1) for a given uq, to label the zeros according to the convention that the 
hrst positive zero is labelled by n = 1. Thus the upper half of the critical line will have the 
zeros labelled by n = 1,2 ,... corresponding to positive tn, while the lower half will have 
the negative values labelled by n = 0, —1,.... The integer uq depends on fc, a and y, 
and should be chosen according to each specihc case. In the cases we analyze below uq = 0, 
whereas for the trivial character no = 1. Henceforth we will omit the integer no in the 
equations, since all cases analyzed in the following have no = 0. Nevertheless, the reader 
should bear in mind that for other cases, it may be necessary to replace n —)■ n — no in the 
following equations. 

In summary, there are an inhnite number of zeros on the critical line, i.e. in the form 
Pn = I + where for a given n G Z, the imaginary part is the solution of the equation 


^k,aitn) + lim arg L (I + (5 + itn, x) - 

5-S.0+ 


arg G (x) 


2 4 


(41) 
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C. Asymptotic equation for the n-th zero 


From Stirling’s formula we have the following asymptotic form for t —>■ ±oo: 


= sgn(t) 


■|t|, fk\t\\ 2a-1 ■ 


( 42 ) 


The hrst order approximation of (41), i.e. neglecting 0{l/t) terms, is therefore given by 
\tr 

-iug \ — - I -r - inn argiv ^ ^ X; 

(43) 


27r 


log lim argL (I+ (5 + iz/„|t„|,x) 

\ 27re / TT 5-s>o+ ^ 

I/n - 4 - 2 a(l + Un) 


- ^ arg G (x) = n + 


where z/„ = 1 if n > 0 and = —1 if n < 0. For n > 0 we have tn = \tn\ and for n < 0 we 
have tn = -\tn\- 

D. Counting formulas 


Let us define iVg'"(T, x) as the number of zeros on the critical line with 0 < Q'(p) < T and 
Nq{T,x) as the number of zeros with —T < 3(p) < 0. As explained before, A"q^(T, y) ^ 
Nq (T, x) if the characters are complex numbers, since the zeros are not symmetrically 
distributed between the upper and lower half of the critical line. 


The counting formula A"q''(T, y) is obtained from equation (41) by replacing —)■ T and 


n —)■ + 4, therefore 


No{T, x) = -^k,a{T) + - argL (4 + iT, x) - ^ argG (x) + y- ( 44 ) 

TT TT 27r 4 


The passage from (41) to (44) is justihed under the assumptions already discussed in con¬ 


nection with (14) and (21), i.e. assuming that (41) has a unique solution for every n. As 


explained above for the Riemann C, case, this is equivalent to assume the existence of the 
hm 5 _,.o+ argL (^ + 5 + R, x) for every t. Analogously, the counting formula on the lower half 
line is given by 

No (T, x) = -^4fc,a(T) - - argL (4 - iT, x) + 7^ argG(x) - y- (45) 

TT TT zn 4 

Note that in (44) and (45) T is positive. Both cases are plotted in Figure |^for the character 


X 7,2 shown in (81). One can notice that they are precisely staircase functions, jumping by 
one at each zero. Note also that the functions are not symmetric about the origin. 
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FIG. 2: Exact counting formulae (44) and (45). Note that they are not symmetric with respect to 


the origin, since the L-zeros for complex y are not complex conjugates. We used y = X 7,2 shown 


in (81). 


From (42) we also have the hrst order approximation for T —)■ co, 


T f kT\ 1 

^ + - argL (i + ^T, x) - 


arg G 1 a 
8 ^ 2 ‘ 


Analogously, for the lower half line we have 


^T- /m \ T f kT\ 1 ,, m \ G 1 

^ - -argL(, -.r,x) + ^ - 


(46) 


(47) 


As in (41) again we are omitting uq since in the cases below uq = 0, but for other cases one 
may need to include ±no on the right hand side of respectively. 

It is known that the number of zeros on the entire critical strip up to height T, i.e. in 
the region {0 < a < 1, 0 < f < T}, is given by [20] 


N^iT.x) = -'&Ka (T) + -argL (I +iT,x) - -argL (|,x) . 

71 71 71 


(48) 


From Stirling’s approximation and noticing that 2a — 1 = —x(—1), then for T —)■ cx) we 
obtain the asymptotic approximation 


(49) 


T f kT\ 1 

^ j +)(^ 0 + 4 

-largL(l,x)-hW + o(l/r). 

TT y 


Both formulas (48) and (49) are exactly the same as (44) and (46), respectively. This 


can be seen as follows. From (30) we conclude that ^ is real on the critical line. Thus 
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arg,^ (|)=0 = —I argG (x) + argL (^, x) + x- Then, replacing argG (x) in ( [4l| ) we obtain 
^k,a {tn) + ^hm argL (^ + h + itn,x) - argL (|,x) = (n - |) vr. (50) 


Replacing tn ^ T and n —)■ Nq + | in (50) we have precisely the expression (48), and also 


(49) for T —)■ cx). Then we conclnde that A^'q^(T, x) = A^+(T, x) exactly. From (31) we see 
that negative zeros for character x correspond to positive zeros for character x- Thus for 
—T < ‘^{p) < 0 the counting on the critical strip also coincides with the counting on the 
critical line, since Nq{T,x) = Nq{T,x) and N~{T,x) = N~^{T,x)- Therefore, the number 
of zeros on the entire critical strip is the same as the number of zeros on the critical line 


obtained as solutions of (41), under the assumption that (41) has a unique solution for every 


n. This is equivalent to stating that lim 5 _,.o+ argL + 6 + it,x) exists for every t. This 
will be further exemplified in Section |VI[ 


IV. ZEROS OF T-FUNCTIONS BASED ON MODULAR FORMS 


Let us generalize the previous results to L-functions based on level one modular forms. 
We first recall some basic definitions and properties. The modular group can be represented 
by the set of 2 x 2 integer matrices 


SL2 {Z) = \A = 


a b 
c d 


a, b,c,d E Z, det A = 1 


(51) 


provided each matrix A is identified with —A, i.e. ±A are regarded as the same transforma¬ 
tion. Thus for T in the upper half complex plane, it transforms as r h- At = under the 
action of the modular group. A modular form / of weight A; is a function that is analytic in 

the upper half complex plane which satisfies the functional relation 

' ar 


f 


* ' = {cT + df f{T). 


(52) 


CT + d^ 

If the above equation is satisfied for all of SL 2 (Z), then / is referred to as being of level 

one. It is possible to define higher level modular forms which satisfy the above equation 

for a subgroup of SL 2 (Z). Since our results are easily generalized to the higher level case, 

henceforth we will only consider level one forms. 

1 


For the SL 2 (Z) element , the above equation (52) implies the periodicity /(r) = 

/(r -h 1), thus it has a Fourier series 


/M = (f. 

n=0 


q = 


(53) 
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If a/(0) = 0 then / is called a cusp form. 

From the Fourier coefficients, one can dehne the Dirichlet series 


L,(s) = 


oo , . 

Q/ jn) 


n=l 




( 54 ) 


The functional equation for Lf (s) relates it to Lf (k — s), so that the critical line is lR(s) = 
where fc > 4 is an even integer. One can always shift the critical line to ^ by replacing 
a/(n) —>■ however we will not do this here. Let us dehne 


Af(s) = {2tt) " F(s) Lf{s). 


(55) 


Then the functional equation is given by 


A/(s) = (-l)‘'"A/(fc-s). 


(56) 


There are only two cases to consider since | can be an even or an odd integer. As in (29) 


we can absorb the extra minus sign factor for the odd case. Thus we dehne if{s) = Af{s) 
for I even, and then ^/(s) = — s). For | odd we dehne ^/(s) = implying 

^f{s) = ^f{k — s). Representing ^f{s) = |'C/| where s = a + it, we follow exactly the same 


steps as in the previous sections. From the particular solution (37) we conclude that there 


are inhnite zeros on the critical line = | determined by lim 5 _,.o+ 0 (| + 5, t) = (n — vr. 
Therefore 
equation 


Therefore, these zeros are given in the form Pn = 2 + where is the solution of the 


^kitn) + ^lim argL/ (I + 5 + R„) = - 


1 + (-1)^/" ^ 


71 


where n = 1,2 ,... and we have dehned 


'dfc(t) = argF (| + it) - t log27r. 


(57) 


(58) 


This implies that the number of solutions of (57) with 0 < t < T is given by 

11 1 — 1—li^A 

No (T) = -MT) + - argL; (I + tT) - 

TT 71 4 

In the limit of large neglecting terms of 0{l/t), the equation ( [^ becomes 

p + (_l)fc/2\ 


(59) 


tn log argT/ (I + 5 + itn) = 


- n. 


(60) 
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V. APPROXIMATE ZEROS IN TERMS OF THE LAMBERT IT-FUNCTION 


A. Explicit formula 

We now show that it is possible to obtain an approximate solution to the previous tran¬ 
scendental equations with an explicit formula. In this approximation, there is indeed a 
unique solution to the equation for every n. Let us introduce the Lambert W-function m, 
which is dehned for any complex number ^ through the equation 


W{z)e 


W(z) ^ 


(61) 


The multi-valued W-function cannot be expressed in terms of other known elementary func¬ 
tions. If we restrict attention to real-valued W(x) there are two branches. The principal 
branch occurs when W{x) > —1 and is denoted by ITo, or simply W for short, and its 
domain is a; > —e~^. The secondary branch, denoted by IT_i, satishes W-i{x) < —1 for 
—e~^ < X < 0. Since we are interested only in positive real-valued solutions, we just need 
the principal branch where W is single-valued. 


Let us start with the zeros of the ^-function, described by equation (13). Consider its 


leading order approximation, or equivalently its average since (arg(^ -|- iy^) = 0. Then we 
have the transcendental equation 


— log! - 
2 ti ^ I 27re 


11 


= n — 


(62) 


Through the transformation tn = 271 {n — this equation can be written as Xn&^" = 


e ^ (n — y). Comparing with (61) we thus we obtain 

27r {n - f) 


W [e-‘ (n - a)] 


(63) 


where n = 1,2,. 


Although the inversion from (62) to (63) is rather simple, it is very convenient since it 


is indeed an explicit formula depending only on n, and W is included in most numerical 
packages. It gives an approximate solution for the ordinates of the Riemann zeros in closed 


form. The values computed from (63) are much closer to the Riemann zeros than Gram 


points, and one does not have to deal with violations of Gram’s law; see Remark 
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Analogously, for Diriclilet L-functions, after neglecting the argL term, the equation (43) 
yields a transcendental equation which can be written as through the 

transformation \tn\ = where 


^ ^ ^ / , argG(x)^^ 1 - 4z/„ - 2a + 1) 

An [x) = i^n[n+ ——— 1 +-^-. 


Thus the approximate solution is explicitly given by 

27TUnAn (x) 


in — 


(64) 


(65) 


HqA^e-M„(x)] 

where n = 0,±1,±2,.... In the above formula n = 1,2,... correspond to positive f, 
solutions, while n = 0, — 1,... correspond to negative solutions. Contrary to the C 
function, in general, the zeros are not conjugate related along the critical line. 


In the same way, ignoring the small argLj term in (60), the approximate solution for the 
imaginary part of the zeros of L-functions based on level one modular forms is given by 


tn 


A„,7r 


IT[(2e)-iA„]’ 


A„, = n- 




n = 1,2,... . 


( 66 ) 


B. Further remarks 


Let us focus on the approximation (63) regarding zeros of the ^-function. Obviously the 


same arguments apply to the zeros of the other classes of functions based on formulas (65) 


and (66) 


Remark 6. The estimates given by (63) can be calculated for arbitrarily large n, since W 


is a standard elementary function. Of course the tn are not as accurate as the solutions 


tn including the arg(^ term, as we will see in Section VII Nevertheless, it is indeed a good 


estimate, especially if one considers very high zeros where traditional methods have not 


previously estimated such high values. For instance, formula (63) can easily estimate the 
zeros shown in Table (Appendix |B 1[ ), and much higher if desirable. The numbers in this 
table are accurate approximations to the n-th zero to the number of digits shown, which 
is approximately the number of digits in the integer part. For instance, the approximation 
to the 10^°° zero is correct to 100 digits. With Mathematica we easily calculated the hrst 


million digits of the 10 


10 ® 


zero. 


Remark 7. Using the asymptotic behaviour W{x) ~ logx for large x, the n-th zero is 
approximately given by tn ~ 27rn/logn, as already known |T0]. The distance between 
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consecutive ordinates is then approximately equal to tn+i — tn ^ 27r/logn, which tends to 
zero when n —)■ cx). 


Remark 8. The solutions (63) are reminiscent of the so-called Gram points which are 


solutions to k){gn) = mi where d is given by (19). Gram’s law is the tendency for Riemann 
zeros to lie between consecutive Gram points, but it is known to fail for about | of all Gram 
intervals. Our are intrinsically different from Gram points. It is an approximate solution 
for the ordinate of the zero itself. In particular, the Gram point go = 17.8455 is the closest 
to the hrst Riemann zero, whereas ti = 14.52 is already much closer to the true zero which is 
ti 14.1347. The traditional method to compute the zeros is based on the Riemann-Siegel 
formula C (^ + = ^{t) [cos'd(f) — zsin'd(f)], and the empirical observation that the real 

part of this equation is almost always positive, except when Gram’s law fails, and Z{t) has 
the opposite sign of sin'd. Since Z(t) and C. + it'j have the same zeros, one looks for 
the zeros of Z{t) between two Gram points, as long as Gram’s law holds {—1)^Z (gn) > 0. 


To verify the RH numerically, the counting formula (22) must also be used to assure that 
the number of zeros on the critical line coincide with the number of zeros on the strip. 
The detailed procedure is throughly explained in mm- Based on this method, amazingly 
accurate solutions and high zeros on the critical line were computed [4HII451 - I5U] . Nevertheless, 


our proposal is fundamentally different. We claim that (20) is the equation that determines 
the Riemann zeros on the critical line. Then, one just needs to hnd its solution for a given 
n. We will compute the Riemann zeros in this way in the next section, just by solving the 


equation (20) numerically, starting from the approximation given by the explicit formula 


(63), without using Gram points nor the Riemann-Siegel Z function. Let us emphasize that 
our goal is not to provide a more efficient algorithm to compute the zeros [IH] , although the 


method described here may very well be, but to justify the validity of equation (20). 


VI. A COUNTEREXAMPLE: THE DAVENPORT-HEILBRONN FUNCTION 

In this section we consider a function that is known to violate the RH, and this serves 
to sharpen our understanding of our previous analysis. In this example, one can clearly see 
how the corresponding transcendental equation does not have a unique solution for every n. 
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The Davenport-Heilbronn function is defined by 


P(s) ^ T 


in) T , . (^ + i^) r r - \ 

. - L (s, X 5 , 2 ) + ^-5- L (s, X'5_2) 


(67) 


with 


K = 


Vl0-2v/5-2 


v/5- 1 


Above the Dirichlet character is the following: 


n 

12 3 

4 5 

X5,2(«) 

1 i —i 

-1 0 


( 68 ) 


(69) 


where X 5 , 2 (—1) = —1 thus a = 1. The function (67) satishes the functional equation 


^( s ) =^(1 - s ), 


m 


7r\ -^/2 


1 + s 


V(s) 


(70) 


The function (67) has almost all the same properties of (, such as a functional equation, 
except that it has no Euler Product Formula. It is well known that it has zeros in the region 
3? (s) > 1, which is essentially a consequence that it has no Euler product. It also has zeros 
in the critical strip 0 < 3fJ(s) < 1, where infinitely many of them lie on the critical line 
3?(s) = however, it also has zeros off of the critical line, thus violating the RH. For a 
detailed study of this function and numerical computation of its zeros see [ST]. 

Repeating the analysis of the previous sections for zeros on the critical line, we obtain 
the following transcendental equation: 


-'d 5 ,i(^n) + - lim argT> (I + 5 + R„) + ^ = 

TT TT 5-S.0+ ^ ^ ^2 


n 


where is dehned in (38). The approximate solution is explicitly given by 

27r (n - I) 


tn. 


W [5 


e~^(n 


( 71 ) 


(72) 


for n = 1,2,.... From (71) and (72) it is possible to compute zeros on the critical line. 


Moreover, zeros off of the critical line satisfy the general solution (10) m- This shows that 


B = 0, with B defined in (35), captures all the zeros. 

Since ( [7l| only captures zeros at a = ^, what happens if there are zeros off of the critical 
line? Consider a simple zero denoted by p, = <J, + it, where 0 < cr, < 1 and cr, 7 ^ |. Due 


to the functional equation (70) there is also a zero at 1 — p, = 1 — a, + it,. Let 


Sv{t) = - lim engV (l + S + it) 

TT (5-S.0+ ^ 


( 73 ) 
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FIG. 3: (a) Left hand side of ( [tT] ) against t. Note the discontinuity at the point t, ^ 85.6993 corre¬ 
sponding to n = 44 and n = 45, where ( pTlj ) has no solution, (b) We plot 1? [V{l/2 + it)] (blue line) 
and O' \p{\l2 -|- it)] (red line) against t. Observe that JR [iD(l/2 -|- it)] < 0 when 3= [V{l/2 -|- it)] —> 0, 
for t —)> t,, signaling the change of branch of argD(l/2 -|- it). 



From its role in the counting formula over the entire critical strip, one knows that when t 
varies across t, then Suit) must jump by two, i.e. we must have ASuit,) = Suit, + e) — 
Suit, — e) =2. This implies that Suit) changes branch around t, in such a way that it 


cannot be smoothed out by the 5 —)■ O’*" limit. In other words, the limit (73) does not exist 


close to t,. Therefore, (71) will not have a solution around t, for a given n. If instead of a 
simple zero we have a zero with multiplicity m > 2, then ASuit,) = 2m, changing branch 
even more drastically. The same situation also happens if there are zeros with multiplicity 
m > 2 on the critical line, where we would have ASu = m. 

In the case of the function ( [6^ , the hrst zero off of the critical line occurs at a, ~ 0.8085 
and t, 85.6993. In Figure we plot the left hand side of (|7l| against t, and one can 


clearly see the above mentioned situation, namely that (71) is not dehned at t, and there 
is no solution for n = 44 and n = 45. The change of branch close to t, can be seen from 
Figure Therefore, denoting iVo(T) the number of solutions of ( [7l| up to height T, we 
clearly have N^iT) < iV(T), where iV(T) is the number of zeros in the entire critical strip. 
For a more detailed illustration of these facts we refer the reader to 


For simple zeros on the critical line the limit (73) exists, but for zeros off of the critical 


line, it does not, since Suit) has to jump at least by two and the change of branch does not 
allow us to smooth the function. 






















(b) 


FIG. 4: Comparison of the prediction of (63) (blue line) and (13) (red dots). We are plotting 
against n. (a) n G [1,...,400]. Note how the solutions are close at first sight, (b) If we focus 


on a small range we can see how the solutions of (13) oscillate around the line (63) due to the 
fluctuating term argC- Here n G [99984,..., 10^]. 

VII. NUMERICAL ANALYSIS: C-FUNCTION 

A. The importance of arg ^ 


Instead of solving the exact equation (20) we will initially consider its first order approxi¬ 


mation, which is equation (13). As we will see, this approximation already yields surprisingly 


accurate values for the Riemann zeros. 


Let us hrst consider how the approximate solution given by (63) is modihed by the 


presence of the arg(C term in (13). Numerically, we compute arg(C taking its principal value. 


The fact that we get very accurate zeros up to the billionth zero implies that up to this t, 
arg ( near a zero is always on the principal branch. As already discussed in Remark [TJ the 
function argC Q + it) oscillates around its average, which is zero, as shown in Figure [^. 
At a Riemann zero it can be dehned by the limit (|^ which is generally not zero. The arg ( 
term plays an important role and indeed improves the estimate of the n-th zero. This can be 
seen in Figure where we compare the estimate given by (63) with the numerical solutions 
of (0. 

We can apply a root hnder method in an appropriate interval, centered around the ap¬ 


proximate solution tn given by formula (63). Some of the solutions obtained in this way 
are 


presented in Table (Appendix B 1) and are accurate up to the number of decimal 
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places shown. We used only Matheniatica or some very simple algorithms to perform these 
numerical computations, taken from standard open source numerical libraries. 


Although the equation (13) was derived for large n, it is surprisingly accurate even for 


the lower zeros, as shown in Table [nT| (Appendix |B 1| ). It is actually easier to solve for low 
zeros since argC is better behaved. These numbers are correct up to the number of digits 
shown, and the precision was improved simply by decreasing the error tolerance. 


B. GUE statistics 


The link between the Riemann zeros and random matrix theory started with the pair 
correlation of zeros, proposed by Montgomery [20] , and the observation of Dyson [56] that it 
is the same as the 2-point correlation function predicted by the Gaussian Unitary Ensemble 
(GUE) for large random matrices. 


The main purpose of this section is to test whether our approximation (13) to the zeros is 


accurate enough to reveal this statistics. Whereas formula (63) is a valid estimate, it is not 


sufficiently accurate to reproduce the GUE statistics, since it does not have the oscillatory 


argC term. On the other hand, the solutions to equation (13) are accurate enough, which 
again indicates the importance of arg^. 


Montgomery’s pair correlation conjecture can be stated as follows: 

1 


N{T) 


0<t,t'<T 


~ /Tub 

J a 


■Jl2y2 


(74) 


where ^ log (D) (« - ('), a < a < j3, N{T) ~ £log(£) aceording to 


and the statement is valid in the limit T —)■ oo. The right hand side of (74) is the 2- 


point GUE correlation function. The average spacing between consecutive zeros is given by 
^ ~ 27r/log (^) —)■ 0 as T —)■ cx). This can also be seen from (63) for very large n, i.e. 
tn+i ——)■ 0 as n —)■ CX). Thus the distance d{t,t') between zeros on the left hand side of 


(74) is a normalized distance. 


While (74) can be applied if we start from the hrst zero on the critical line, it is unable 


to provide a test if we are centered around a given high zero on the line. To deal with such 
a situation, Odlyzko [50] proposed a stronger version of Montgomery’s conjecture by taking 
into account the large density of zeros higher on the line. This is done by replacing the 
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M = 1, N = 10 ® 


M = 10® - 10®, Af = 10® 




FIG. 5: The solid line represents the right hand side of (76) and the dots represent its left hand 

, 3) and 


side, computed from equation ( [l^ . The parameters are /3 = a + 0.05, a = (0, 0.05, 
the x-axis is given by x = ^ (a + /3). (a) We use the hrst 10® zeros, (b) The same parameters but 
using zeros in the middle of the critical line; M = 10® — 10® and N = 10®. 


normalized distance in (74) by a sum of normalized distances over consecutive zeros in the 
form 

1 / y- \ 

(75) 


dn — 2^ log f 27 ,- ) ■ 


Thus (74) is replaced by 

1 


E' 


f/3 


1 = 


du I 1 


sm TTM 


■J^2y2 


(76) 


(N - M)(l}-a) ^ ■ B-a 

where M is the label of a given zero on the line and N > M. In this sum it is also assumed 


that n > m, and we included the correct normalization on both sides. The conjecture (76) 
is already well supported by extensive numerical analysis 


Odlyzko’s conjecture (76) is a very strong constraint on the statistics of the zeros. Thus we 


submit the numerical solutions of equation (13) to this test. In Figure]^ we can see the result 
for M = 1 and N = 10®, with a ranging from 0 ... 3 in steps of e = 0.05, and jd = a+e for each 
value of a, i.e. a = (0.00, 0.05, 0.10,..., 3.00) and fd = (0.05, 0.10, ..., 3.05). We compute 


the left hand side of (76) for each pair {a, fd) and plot the result against x = \ + (d). 
In Figure we do the same thing but with M = 10® — 10® and N = 10®. Clearly, the 
numerical solutions of (13) reproduce the CUE statistics. In fact, Figure]^ is identical to 
the one in [50] . The last zeros in these ranges are shown in Table |IV| (Appendix |B 1[) . 
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C. Prime number counting function 


In this section we explore whether our approximations to the Riemann zeros are accurate 
enough to reconstruct the prime number counting function. As usual, let 7r(x) denote the 
number of primes less than x. Riemann obtained an explicit expression for 7r(x) in terms of 
the non-trivial zeros of ((s). There are simpler but equivalent versions of the main result, 
based on the function ip{x) below. However, let us present the main formula for 7i{x) itself 
since it is historically more important. 

The function 7i{x) is related to another number-theoretic function J{x), dehned as 


J{x) 


E 

2<n<x 


A(n) 

logn 


(77) 


where A(n), the von Mangoldt function, is dehned as A(n) = logp if n = for some prime 
p and an integer m > 1, and A(n) = 0 otherwise. The two functions 7r{x) and J{x) are 
related by Mbbius inversion: 

’^(^1 = E (78) 

n>l 

Here p{n) is the Mobius function dehned as follows. p(n) = 0 if n has one or more repeated 
prime factors, /i(n) = 1 if n = 1 and /i(n) = (—1)^ if n is a product of k distinct primes. 
The above expression is actually a hnite sum, since for large enough n, x^^^ < 2 and J = 0. 

The main result of Riemann is a formula for J [x) , expressed as an inhnite sum over zeros 
p of the C(-s) function: 

J(x) = U(x)~^U(xn + l pi7VT)-'°8-2. (79) 

where Li(a;) = f^dt/log t is the log-integral function^. The above sum is real because the 
p’s come in conjugate pairs. If there are no zeros on the line = 1, then the dominant 
term is the hrst one in the above equation, J(x) ~ Li(a;), and this was used to prove the 
prime number theorem by Hadamard and de la Vallee Poussin. 

The function ip{x) has the simpler form 

i’ix) = = ~ ^ log (80) 

n<x P ^ \ ^ / 


^ Some care must be taken in numerically evaluating Li(a;^) since Li has a branch point. It is more properly 
defined as Ei(plogx) where Ei(z) = — dt e~^/t is the exponential integral function. 
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FIG. 6: The prime number counting function 'k{x) with the first 50 Riemann zeros, (a) Zeros 


approximated by the formula (63). (b) Zeros obtained from numerical solutions to the equation 


(13). 


In this formulation the prime number theorem is equivalent to 'ip{x) ~ x. 

In Figurewe plot 7i{x) from equations (78) and (79), computed with the hrst 50 zeros 
in the approximation Pn = | + itn given by (63). Figure & shows the same plot with zeros 
obtained from the numerical solutions of equation (13). Although with the approximation 
tn the curve is trying to follow the steps in ^{x), once again, one clearly sees the importance 
of the arg ( term. 


D. Solutions to the exact equation 


In the previous sections we have computed numerical solutions of (13) showing that. 


actually, this hrst order approximation to (20) is very good and already captures some 


interesting properties of the Riemann zeros, such as the GUE statistics and the ability to 


reproduce 7r{x). Nevertheless, by simply solving (20) it is possible to obtain values for the 
zeros as accurately as desirable. The numerical procedure is performed as follows: 


1. We apply a root hnder method on (20) looking for the solution in a region centered 
around the number tn provided by (63), with a not so small 6, for instance 6 ~ 10“®. 


2. We solve (20) again but now centered around the solution obtained in step above, 
and we decrease 6, for instance 6 ~ 10“®. 


3. We repeat the procedure in step above, decreasing 6 again. 
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4. Through successive iterations, and decreasing S each time, it is possible to obtain 
solutions as accurate as desirable. In carrying this out, it is important to not allow 6 
to be exactly zero. 


An actual implementation of the above procedure in Mathematica is shown in Appendix 
which we have included mainly to show its simplicity. The first few zeros computed in 
this way are shown in Table [V| (Appendix |B 1[ ). Through successive iterations it is possible 
achieve even much higher accuracy than shown in Table |Vj 

It is known that the hrst zero where Gram’s law fails is for n = 126. Applying the same 


method, like for any other n, the solution of (20) starting with the approximation (63) does 
not present any difficulty. We easily found the following number: 


fi26 = 279.229250927745189228409880451955359283492637405561293594727 


Just to illustrate, and to convince the reader, how the solutions of (20) can be made arbi 


trarily precise, we compute the zero n = 1000 accurate up to 500 decimal places, also using 
the same simple approach'^: 


Gooo = 1419.42248094599568646598903807991681923210060106416601630 
46908146846086764175930104179113432911792099874809842 
32260560118741397447952650637067250834288983151845447 
68825259311594423942519548468770816394625633238145779 
15284185593431511879329057764279980127360524094461173 
37041818962494747459675690479839876840142804973590017 
35474131911629348658946395454231320810569901980719391 
75430299848814901931936718231264204272763589114878483 
29996467356160858436515425171824179566414953524432921 
93649483857772253460088 


Furthermore, one can substitute known precise Riemann zeros into (20) and can check that 


the equation is identically satished. These results corroborate that (20) is an exact equation 
for the Riemann zeros. 


Computing this number to 500 digit accuracy took a few minutes on a standard personal laptop computer 
using Mathematica. It only takes a few seconds to obtain 100 digit accuracy. 
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VIII. NUMERICAL ANALYSIS: L-FUNCTIONS 


We perform exactly the same numerical procedure as described in the previous Sec¬ 


tion VIID, but now with equation (41) and (65) for Dirichlet L-functions, or with (57) and 


( 66 ) for L-functions based on level one modular forms. 


A. Dirichlet L-functions 


We will illustrate our formulas with the primitive characters X 7,2 and Xi,^ since they 
possess the full generality of a = 0 and a = 1 and complex components. There are actually 
99 ( 7 ) = 6 distinct characters mod 7. 

Example X 7 , 2 - Consider k = 7 and j = 2, i.e. we are computing the Dirichlet character 
X 7 2 (n). For this case a = 1. Then we have the following components: 


(81) 


n 

1 2 3 4 5 

6 7 

XtA^) 

g27ri/3 g7ri/3 g—27ri/3 ^—nil3 

1 

0 


The hrst few zeros, positive and negative, obtained by solving (41) are shown in Table VI 


m 


Appendix B 2 The solutions shown are easily obtained with 50 decimal places of accuracy. 


Example X 7 , 3 - Consider k = 7 and j = 3, such that a = 0. In this case the components of 
X 7 , 3 (’^) S'l’e the following: 


n 


X7M 


6 7 


(82) 


X g-2m/3 g27ri/3 g27ri/3 ^-2m/3 ^ g 

The hrst few solutions of (41) are shown in Table VII in Appendix |B 2| and are accurate up 
to 50 decimal places. As previously stated, the solutions to equation ( [4l| can be calculated 
to any desired level of accuracy. For instance, we can easily compute the following number 
for n = 1000 , accurate to 100 decimal places: 

= 1037.56371706920654296560046127698168717112749601359549 
01734503731679747841764715443496546207885576444206 

We also have been able to solve the equation for high zeros to high accuracy, up to the 
millionth zero, some of which are listed in Table VIII in Appendix |B 2[ and were previously 
unknown. 
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B. Modular i-function based on Ramanujan r 


Here we will consider an example of a modnlar form of weight k = 12. The simplest 
example is based on the Dedekind 77 -fnnction 

00 

rt{T) = n(l - «”)• 9 = 

n=\ 


Up to a simple factor, rj is the inverse of the chiral partition fnnction of the free boson 
conformal held theory [5l], where r is the modnlar parameter of the torns. The modnlar 
discriminant 

00 

A(T) = ,,(ry‘ = ^T(n),” (84) 

n=l 

is a weight k = 12 modnlar form. It is closely related to the inverse of the partition fnnction 
of the bosonic string in 26 dimensions, where 24 is the nnmber of light-cone degrees of 
freedom in 26 spacetime dimensions [55] • The Fonrier coefficients r(n) correspond to the 
Ramannjan r-fnnction, and the hrst few are 


n 

12 3 

4 

5 

6 

7 

8 

r(n) 

1 -24 252 

-1472 

4830 

-6048 

-16744 

84480 


We then dehne the Dirichlet series 


n=\ 




Applying (57) the zeros are p„ 


6 -f- where tn satishes the eqnation 


(85) 


( 86 ) 


^i 2 (tn) + ^lim arg La(6 + 5 + itn) — (n - 4] tt. 


The counting formula (59) and its asymptotic approximation are 

7Vo(T) = Uu{T) + - argLA(6 + tT) 

71 71 

T, / T \ 1 , ,11 

- - log U— + -argLA(6 + «T) + — . 
71 V 27re J 71 4 


(87) 

( 88 ) 
(89) 


A plot of ( 88 ) is shown in Fignrej^ and we can see that it is a perfect staircase fnnction 
The approximate solntion ( | 6 ^ now has the form 

(" - f) 


trr. 


TT 


W [(2 


e)~^ (n 


¥)] 


for n = 2,3, 


(90) 
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FIG. 7: Exact counting formula (88) based on the Ramanujan r-function. 


Note that the above equation is valid for u > 1 since W{x) is not dehned for x < —1/e. 


We follow exactly the same numerical procedure, previously discussed in Section |VII D 
and implemented in Appendix]^ to solve the equation (87) starting with the approximation 


provided by (90). Some of these solutions are shown in Table IX in Appendix B3 and are 


accurate to 50 decimal places. 


IX. CONCLUDING REMARKS 


In this paper we considered non-trivial zeros of the Riemann (/-function, Dirichlet L- 
functions and L-functions based on level one modular forms. The same approach was applied 
to all these classes of functions, showing that there are an inhnite number of zeros on the 
critical line in one-to-one correspondence with the zeros of the cosine function ( pT| , leading 
to a transcendental equation satished by the ordinate of the n-th zero. More specihcally, for 


the Riemann (/-function these zeros are solutions to (20), for Dirichlet L-functions we have 


(41), and for L-functions based on level one modular forms the ordinates of the zeros must 


satisfy (57). It is important to stress that these equations were derived on the critical line, 


without assuming the RH. 

The implication of our work for the GRH can be summarized as follows. If the corre¬ 
sponding transcendental equation has a unique solution for every n, the validity of the GRH 
would follow. The explanation is very simple. Suppose the transcendental equation indeed 
has a unique solution for every n. Then the zeros obtained from its solutions on the critical 
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line can be counted, since they are enumerated by the integer n, yielding the counting func¬ 
tion Nq{T). The number of solutions saturate the counting formula over the entire critical 
strip, namely Nq{T) = N{T), where N{T) counts zeros on the entire critical strip and has 
been known for a long time. Thus the equation captures all the non-trivial zeros. 


As previously discussed, and explicitly illustrated in Section |VTJ the existence of solutions 
depends on whether the 5 —)■ O’*" limit of the argument of the corresponding L-function is well- 
dehned for every ordinate t. The validity of this limit was our only assumption throughout 
the paper. We also argued that if there is indeed a unique solution of the transcendental 
equation for every n, then all non-trivial zeros are simple. If there are zeros off of the 
critical line, or zeros with multiplicity m > 2 on the critical line, the equation will fail to 
capture all the zeros on the critical strip and then No{T) < N{T). Does this means that the 
GRH is false if the transcendental equation does not always have a unique solution? Not 
necessarily, since all the zeros can still be on the critical line but not all of them are simple. 
This suggests that the GRH and the simplicity of all non-trivial zeros is equivalent to the 
statement that the transcendental equation has a unique solution for every n. We have not 
proven that there is a unique solution to the transcendental equation. An attempt to justify 
more carefully this 5 —)■ O’*" limit is in our preliminary work [58], where we claim that the 
Euler Product Formula is still valid in the region | < 3fJ(s) < 1 in a statistical manner. 

We also showed that it is possible to obtain an explicit formula as an approximation for 


the ordinates of the zeros in terms of the Lambert W-function; equation (63) for the C- 


function, (65) for Dirichlet L-functions and (66) for L-functions based on level one modular 
forms. This approximation is very convenient, allowing us to actually compute accurate 
zeros without relying on Gram points, nor dealing with violations of Gram’s law. 

We have also provided compelling numerical evidence for the validity of these transcen¬ 
dental equations satished by the n-th zero. For the (^-function, the leading order asymptotic 


approximation (13) proved to be accurate enough to reveal the interesting features of the 
Riemann zeros, like the GUE statistics and the reconstruction of the prime number counting 


function 7r(x). It turns out the exact equation (20) is much more stable and easy to solve 
numerically, it is thus able to provide numerical results as accurate as is desired. We have 


also provided accurate numerical solutions for Dirichlet L-functions using (41) and for the 
particular example of the modular L-function based on the Ramanujan r-function, through 


(87). The numerical approach employed here constitutes a novel and simple method to 
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compute non-trivial zeros of L-functions. 
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Appendix A: Mathematica implementation 


Here we provide the short Mathematica code used to compute the zeros from the transcen¬ 
dental equations. We will consider Dirichlet L-functions, since it involves more ingredients, 
like the modulus k, the order a and the Gauss sum G{t). For the Riemann (^-function 
the procedure below is trivially adapted as a special case, as well as for the Ramanujan 
r-function of Section IVHl B[ 


The function (38) is implemented as follows: 


RSTheta[t_, a_, k_] := Im[LogGamma[l/4-|-a/2-|-I*t/2]] — t/2*Log[Pi/k] 

For the transcendental equation ( 

41 

) we have 


ExactEq[n_, t_, s_, a_, k_, j_, G_, n0_] := 

(RSTheta[t, a, k] -I- Arg[DirichletL[k, j, l/2-|-\delta-|-I*t]] — l/2*Arg[G])/Pi -I- a/4 -|- 1/2 — n 
nO 


Above, s denotes 0 < 5 <C 1, a is the order (26), k is the modulus, j specihes the Dirichlet 


character Xk,i (as discussed in Section HI), and G is the Gauss sum (25). Note that we 


also included no, discussed after (40), but we always set uq = 0 for the cases analysed in 


Section VlllA The implementation of the approximate solution (65) is 


Sgn[n_] := Which[n != 0, Sign[n], n == 0, —1] 

A[n_, a_, G_, n0_] := Sgn[n]*(n — nO -I- l/2/Pi*Arg[G]) 

-I- (1 — 4*Sgn[n] — 2*a*(Sgn[n]-|-l))/8 
tApprox[n_, a_, G_, k_, n0_] := 

2*Pi*Sgn[n]*A[n, a, G, nO]/LambertW[k*A[n, a, G, nO]/E] 


One can then obtain the numerical solution of the transcendental equation (41) as follows: 
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FindZero[n_, s_, a_, k_, j_, G_, nO_, tO_, prec_] := 

t /. FindRoot[ExactEq[n, t, s, a, k, j, G, nO], {t, tO}, PrecisionGoal— >prec/2, AccuracyGoal 
— >prec/2, WorkingPrecision— >prec] 


Above, to will be given by the approximate solution (65). The variable prec will be adjusted 


iteratively. Now the procedure described in Section VIID can be implemented as follows: 


DirichletNZero[n_, order., digits., k., j., nO.] := ( 
chi = DirichletCharacter[k, j, —1]; 
a = Which[chi == —1, 1, chi == 1, 0]; 
s = lOA-S); 
prec = 15; 

G = Sum[DirichletCharacter[k, j, l]*Exp[2*Pi*I*l/k], {1, 1, k}]; 
t = N[tApprox[n, a, G, k, nO], 20]; 
absvalue = 1; 

While[absvalue > order, 

t = FindZero[n, s, a, k, j, G, nO, t, prec]; 
Print[NumberForm[t, digits]]; 
s = s/1000; 
prec = prec + 20; 

absvalue = Abs[DirichletL[k, j, 1/2 + I*t]]; 

] 

Print [ScientificForm[absvalue, 5]]; 


Above the variable order controls the accuracy of the solution. For instance, if order = 
it will iterate until |L (|) + it) | ~ 10“®°. The variable digits controls the number of decimal 
places shown in the output. 


Let us compute the zero n = 1, for the character (82), i.e. k = 7 and j = 3. We will 
verify the solution to order = 10“^° and print the results with digits = 22. Thus executing 


DirichletNZero]!, 10^-20), 22, 7, 3, 0] 


the output will be 
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4.35640188194944 

4.356401624736541498075 

4.356401624736284227537 

4.356401624736284227280 

4.1664*10^(-25) 


Note how the decimal digits converge in each iteration. By decreasing order and increasing 
digits it is possible to obtain highly accurate solutions. It is exactly in this way that we 
obtained the tables shown in Appendix Obviously, depending on the height of the critical 
line under consideration, one should adapt the parameters s and prec appropriately. In 
Mathematica we were able to compute solutions up to n ~ 10® for Dirichlet L-functions, 
and up to n ~ 10® for the Riemann (^-function without problems. We were unable to go 
much higher only because Mathematica could not compute the arg L term reliably. To solve 


the transcendental equations (20) and (41) for very high values on the critical line is still a 
challenging numerical problem. Nevertheless, we believe that it can be done through a more 
specialized implementation. 


Appendix B: Numerical results 

In this section we present some of the numerical results obtained by solving the tran¬ 
scendental equations described in this paper. The numerical procedure is described in Ap¬ 
pendix 1^ and should be adapted to each particular class of functions. 


1. Riemann ^-function 


The explicit formula (63) can estimate very high Riemann zeros, yielding results accurate 
up to the decimal point. Some of these results are shown below: 


n tn 

10^2 1.370919909931995308226636 xlO^^ 

10^° 5.741532903784313725642221053588442131126693322343461 x10^® 

(continued) 
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(continued) 

IQi™ 2.806903838428940699031954458382564000845480301628460 

45192360059224930922349073043060335653109252473234 xlO^® 
lO^oo 1.385792222146789340845466805467159190123402451538707 

081832868352483938909689796343076797639408172610028 
651791994879400728026863298840958091288304951600695 
814960962282888090054696215023267048447330585768 xlO^®* 


TABLE I: Numerical results predicted by formula (63), which can easily estimate very high Rie- 
mann zeros. The results are expected to be correct up to the decimal point, i.e. to the number of 
digits in the integer part. The numbers are shown with three digits beyond the integer part. 


In the following table we have numerical solutions to (13), obtained simply by applying 


a root hnder method around the estimate provided by formula (63): 


n 



1 

14.52 

14.134725142 

10 

50.23 

49.773832478 

10^ 

235.99 

236.524229666 

10^ 

1419.52 

1419.422480946 

10-* 

9877.63 

9877.782654006 

10^ 

74920.89 

74920.827498994 

10® 

600269.64 

600269.677012445 

10" 

4992381.11 

4992381.014003179 

10® 

42653549.77 

42653549.760951554 

10^ 

371870204.05 

371870203.837028053 

lOlo 

3293531632.26 

3293531632.397136704 


TABLE IL Numerical solutions to the asymptotic equation (13). All numbers shown are accurate 


up to the 9-th decimal place and agree with |48| . 


Decreasing the error tolerance we can obtain more accurate solutions to the asymptotic 


equation (13), even for the lower zeros, as shown below: 
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n tn 

1 14.13472514173469379045725198356247 

2 21.02203963877155499262847959389690 

3 25.01085758014568876321379099256282 

4 30.42487612585951321031189753058409 

5 32.93506158773918969066236896407490 


TABLE III: Numerical solutions to (13) for the lowest zeros. Although it was derived for high t, it 


provides accurate solutions even for the lower zeros. These numbers are correct up to the decimal 
place shown 


While the previous tables were computed for isolated zeros, to test Odlyzko-Montgomery 


pair correlation conjecture (76) we have to compute systematically a wide range of zeros. 


This is a strong test of equation (13) and the approximation (63), since in principle it could 
have missed some zeros or presented some numerical issues. This was dehnitely not the case. 
We computed all the zeros in the range n = 1... 10^ and also n = 10® — 10^ ... 10®. The 


equation (13) and also the approximation (63) did not miss a single zero. The last numbers 
in these ranges are shown below: 


71 


n 


10^ - 5 

74917.719415828 

10^ - 5 

371870202.244870467 

10^-4 

74918.370580227 

10^-4 

371870202.673284457 

10^-3 

74918.691433454 

10^ - 3 

371870203.177729799 

10^-2 

74919.075161121 

10^ - 2 

371870203.274345928 

lO'^ - 1 

74920.259793259 

10^ - 1 

371870203.802552324 

lO'^ 

74920.827498994 

10^ 

371870203.837028053 

TABLE IV: Last numerical solutions to (13) around n 

= 10^ and n = 10®. 


The previous numerical solutions to (13) were obtained with no iteration, i.e. simply by 
applying the root hnder function once. 


The numerical solutions to the exact equation (20) can yield arbitrarily accurate values. 
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With some very few iterations, as described in Appendix]^ we compnted the hrst few zeros: 


n tn 

1 14.1347251417346937904572519835624702707842571156992431756855 

2 21.0220396387715549926284795938969027773343405249027817546295 

3 25.0108575801456887632137909925628218186595496725579966724965 

4 30.4248761258595132103118975305840913201815600237154401809621 

5 32.9350615877391896906623689640749034888127156035170390092800 


TABLE V: The first few numerical solutions to (20), accurate to 60 digits (58 decimals). 


Dirichlet L-functions 


Below we present some nnmerical solntions to (41), with the Dirichlet character shown 


in (81). We used exactly the procedure described in Appendix [A} For n > 0 we have the 
zeros on the upper half of the critical line, while for n < 0 we have the zeros on the lower 
half of the critical line. 


n 


tri 




10 25.57 25.68439458577475868571703403827676455384372032540097 
9 23.67 24.15466453997877089700472248737944003578203821931614 
8 21.73 21.65252506979642618329545373529843196334089625358303 
7 19.73 19.65122423323359536954110529158230382437142654926200 
6 17.66 17.16141654370607042290552256158565828745960439000612 
5 15.50 15.74686940763941532761353888536874657958310887967059 
4 13.24 13.85454287448149778875634224346689375234567535103602 
3 10.81 9.97989590209139315060581291354262017420478655402522 
2 8.14 8.41361099147117759845752355454727442365106861800819 
1 4.97 5.19811619946654558608428407430395403442607551643259 

(continued) 
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(continued) 


0 

-3.44 

-1 

-7.04 

-2 

-9.85 

-3 

-12.35 

-4 

-14.67 

-5 

-16.86 

-6 

-18.96 

-7 

-20.99 

-8 

-22.95 

-9 

-24.87 


-2.50937455292911971967838452268365746558148671924805 
-7.48493173971596112913314844807905530366284046079242 
-9.89354379409772210349418069925221744973779313289503 
-12.25742488648921665489461478678500208978360618268664 
-14.13507775903777080989456447454654848575048882728616 
-17.71409256153115895322699037454043289926793578042465 
-18.88909760017588073794865307957219593848843485334695 
-20.60481911491253262583427068994945289180639925014034 
-22.66635642792466587252079667063882618974425685038326 
-25.28550752850252321309973718800386160807733038068585 


TABLE VI: Numerical solutions to (41) starting with the approximation (65), for the character 


(81). The solutions are accurate to 50 decimal places. 


The next two tables contains numerical solutions to (41), but with the Dirichlet character 


(82). 


10 25.55 

9 23.65 

8 21.71 

7 19.71 

6 17.64 

5 15.48 

4 13.21 

3 10.79 

2 8.11 
1 4.93 

0 -5.45 
-1 -8.53 


26.16994490801983565967242517629313321888238615283992 

23.20367246134665537826174805893362248072979160004334 

21.31464724410425595182027902594093075251557654412326 

20.03055898508203028994206564551578139558919887432101 

17.61605319887654241030080166645399190430725521508443 

15.93744820468795955688957399890407546316342953223035 

12.53254782268627400807230480038783642378927939761728 

10.73611998749339311587424153504894305046993275660967 

8.78555471449907536558015746317619235911936921514074 

4.35640162473628422727957479051551913297149929441224 

-6.20123004275588129466099054628663166500168462793701 

-7.92743089809203774838798659746549239024181788857305 


(continued) 
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(continued) 


-2 

-11.15 

-3 

-13.55 

-4 

-15.80 

-5 

-17.94 

-6 

-20.00 

-7 

-22.00 

-8 

-23.94 

-9 

-25.83 


-11.01044486207249042239362741094860371668883190429106 

-13.82986789986136757061236809479729216775842888684529 

-16.01372713415040781987211528577709085306698639304444 

-18.04485754217402476822077016067233558476519398664936 

-19.11388571948958246184820859785760690560580302023623 

-22.75640595577430793123629559665860790727892846161121 

-23.95593843516797851393076448042024914372113079309104 

-25.72310440610835748550521669187512401719774475488087 


TABLE VII: Numerical solutions of (41) starting with the approximation (65), for the character 


(82). The solutions are accurate to 50 decimal places. 


n tn 

10^ 1037.61 

10^ 7787.18 

lO*^ 61951.04 

10® 512684.78 


1037.563717069206542965600461276981687171127496013595490 

7787.337916840954922060149425635486826208937584171726906 

61950.779420880674657842482173403370835983852937763461400 

512684.856698029779109684519709321053301710419463624401290 


TABLE VIIL Higher zeros for the Dirichlet character (82). These solutions to (41) are accurate 


to 50 decimal places. 


3. L-function based on Ramanujan r 

Adapting the numerical procedure ofj^for the modular L-function based on the Ramanu¬ 
jan r-function, described in section |VlIIB[ we can obtain the following numerical solutions, 
some of which were previously unknown: 


9.22237939992110252224376719274347813552877062243201 

13.90754986139213440644668132877021949175755235351449 

(continued) 


1 

2 12.46 
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(continued) 


3 

16.27 

4 

19.30 

5 

21.94 

6 

24.35 

7 

26.60 

8 

28.72 

9 

30.74 

10 

32.68 

100 

143.03 

200 

235.55 

300 

318.61 


17.44277697823447331355152513712726271870886652427527 

19.65651314195496100012728175632130280161555091200324 

22.33610363720986727568267445923624619245504695246527 

25.27463654811236535674532419313346311859592673122941 

26.80439115835040303257574923358456474715296800497933 

28.83168262418687544502196191298438972569093668609124 

31.17820949836025906449218889077405585464551198966267 

32.77487538223120744183045567331198999909916163721260 

143.08355526347845507373979776964664120256210342087127 

235.74710143999213667703807130733621035921210614210694 

318.36169446742310747533323741641236307865855919162340 


TABLE IX: Non-trivial zeros of the modular L-function based on the Ramanujan r-function, 


obtained from (87) starting with the approximation (90). 


These solutions are accurate to 50 


decimal places. 
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